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PREFACE 


The  investigation  reported  herein  was  conducted  in  response  to  a 
request  from  the  U.  S.  Army  Engineer  District,  Hew  Orleans  (liOD).  Funds 
were  authorized  by  NOD  on  3 January  1977. 

The  study  was  conducted  by  Dr.  John  J.  Wanstrath  of  the  Coastal 
Branch,  Wave  Dynamics  Division,  Hydraulics  Laboratory,  U.  S.  Army  Engi- 
neer Waterways  Experiment  Station  (WES),  under  the  general  supervision 
of  Dr.  C.  L.  Vincent,  Chief  of  the  Coastal  Branch,  Dr.  R.  W.  Whalin, 
Chief  of  the  Wave  Dynamics  Division,  and  Mr.  H.  B.  Simmons,  Chief  of 
the  Hydraulics  Laboratorj'. 

Director  of  WES  during  the  conduct  of  the  investigation  and  the 
preparation  and  publication  of  this  report  was  COL  John  L.  Cannon,  CE. 
Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  ISITRIC  (Sl)  TO  U.  S.  CUSTOMARY  AKD 
U.  S.  CUSTOMARY  TO  I-STRIC  (Sl)  UNITS  OF  MEASUREMENT 


Units  of  measurement  used  in  this  report  can  he  converted  as  fcllovs; 


Multiply 

By 

To  Obtain 

Metric 

(SI)  to  U.  S. 

Customary 

metres 

3.280839 

feet 

kilometres 

0.5399568 

miles  (u.  S.  nautical) 

metres  per  second 

3.280839 

feet  per  second 

square  centimetres 
per  second 

0.1550 

square  inches  per  second 

millibars 

0.011+50377 

pounds  per  square  inch 

U.  S. 

Customary  to  Metric  (Sl) 

inches 

25.1+ 

millimetres 

feet 

0.301+8 

metres 

miles  (U.  S.  nautical) 

1.852 

kilometres 

miles  (U.  S.  statute) 

1.60931+1+ 

kilometres 

square  feet 

O.0929030I+ 

square  m.etres 

cubic  feet 

0.02831685 

cubic  metres 

pounds  (force)  per 
square  inch 

669I+.757 

pascals 

feet  per  second 

O.30I+8 

metres  per  second 

miles  per  hour 
(U.  S.  statute) 

1.60931+1+ 

kilometres  per  hour 

degrees  (angle) 

0,0171^5329 

radians 

m OPEN-COAST  MATHgl^TICi^-L  STORM  SURGE  MODEL  WITH 
COASTAL  FLOODING  FOR  LOUISIMA 

THEORY  AID  APPLICATION 

PART  I : INTRODUCTION 

1.  Three-dimensional  partial  differential  equations  govern  the 
motion  of  an  infinitesimal  fluid  element.  These  equations  result  from 
hasic  considerations  of  mass  conservation  and  Nevton's  second  law  of 
motion.  The  assumptions  involving  the  incompressihility  and  homogeneity 
of  the  water,  negligihle  vertical  accelerations  of  the  fluid  parcel, 
reasonably  uniform  horizontal  flow  over  the  fluid  depth,  and  others  re- 
sult in  the  classical,  vertically  integrated  long-wave  equations  of 
motion  and  mass  conservation.  Energy  is  supplied  at  the  free  surface, 
in  general,  through  the  action  of  the  wind  and  dissipated  at  the  sea- 
bed through  friction.  The  system  of  equations  is  time-dependent,  two- 
dimensional  in  terms  of  the  horizontal  coordinates,  and  readily  in- 
tegrable  through  numerical  techniques.  These  equations  are  applicable 
to  the  study  of  storm  surge  generation^  on  the  continental  shelf  and 

to  the  degree  that  the  assumptions  are  valid,  nearshore  astronomical 
tide  simulation. 

2.  In  the  past,  numerical  integration  of  the  two-dimensional  as 

well  as  the  three-dimensional  equations  of  motion  have  been  performed 
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using  rectilinear  grids.  Recently,  curvilinear  coordinates  ’ ’ have 
been  employed  in  two-dimensional  models,  with  particular  applications 
to  free  and  forced  long-wave  simiulations  for  large  (hundreds  of  miles*) 
open-coast  stretches  of  the  continental  shelf. 

3.  Open-coast  curvilinear  models  are  considered  superior  to 
rectilinear  m.odels  because  rectilinear  models  must  represent  the  coast 
boundary  as  a series  of  straight-line  segments  connected  at  right  angles. 

* A table  of  factors  for  converting  metric  (3l)  units  of  measurement  to 
U,  S.  customary  units  and  U.  S.  customary  units  to  metric  (Sl)  units 
is  presented  on  page  3. 


Spurious  oscillations  are  injected  into  the  calcxilation  by  the  boundary. 
Furthermore,  the  stair-step  boundary  must  retain  more  vater  over  that 
natural  coastal  configuration  where  water  is  free  to  move  artificially 
unobstructed  and  in  accord  with  the  forces. 

The  use  of  any  of  the  previously  mentioned  open-coast  models 
in  simulating  circulation  and  water  level  conditions  is  made  difficult 
because  of  the  requirement  of  specifying  appropriate  boundary  conditions. 
This  is  especially  the  case  in  the  nearshore  region  where  these  models 
simulate  the  coast  boundary  as  an  infinitely  high,  continuous  wall. 

This  boundary  condition  neglects  overtopping  of  low-lying  land  and  bay 
communication  with  the  open  sea.  These  coastal  processes  have  a con- 
siderable effect  on  the  nearshore  water  levels  and  fluid  velocities. 
Furthermore,  it  is  precisely  this  zone  where  much  environmental  interest 
is  focused  and  the  use  of  numerical  models  should  prove  beneficial. 

5.  This  report  describes  a more  appropriate  open-coast  boundary'’ 
condition  and  presents  results  from  foiir  historical  hurricane  surges 
which  affected  the  Louisiana  coast  from  Atchafalaya  Bay  to  the  Missis- 
sippi River.  The  boundary  conditions  is  termed  the  Finite  Height  Barrier 
Coast  (FHBC)  and  is  incorporated  into  a two-dimensional  model  that 
employs  an  orthogonal  curvilinear  coordinate  system  with  telescoping 
computing  cells.  This  hydrodynamic  model  with  the  FHBC  is  named 
SSURGE  III.  In  a companion  report,^  results  from  other  verification 
studies  of  SSURGE  III  are  presented, 

6.  Three  numerical  programs  are  required  for  staging  the 
production  of  surge  or  tide  computations.  The  first  numerical  program 
determines  the  conformal  mapping  coefficients  for  the  particular  study 
region  of  the  continental  shelf.  The  actual  coastline,  as  seen  on 
national  Ocean  Survey  (NOS)  nautical  charts,  is  smoothed  relative  to 
portraying  small-sca,le  features.  This  continuous  coastline  and  a deep- 
sea  boundary  curve  (perhaps,  following  the  300-  or  600-ft  isobath)  are 
the  input  to  the  conformal  mapping  program.  The  coefficients  are  de- 
termined and  input  to  the  second  numerical  program  which  determines  the 
computing  grid  data  relative  to  a particular  design  by  the  user.  The 
grid  data  are,  in  part,  the  input  to  the  SSURGE  III  program.  Other 


input  includes  the  shelf  bathjr-metry,  FHBC  data,  hydrograph  and  velocity 
output  locations,  and  other  readily  determinable  parameters.  In  actual 
practice,  the  entire  process  is  not  time— consuming  or  expensive. 


PART  II:  COIIFjRMAL  MAPPING 


7.  It  is  desired  to  confomally  map  a spatial  region  of  prototype 
space  of  the  continental  shelf  into  a rectangle  in  a mathematical  image 
plane  in  which  the  coastline  and  deep-sea  boundary  curves  are  specifi- 
cally transformed  into  the  image  plane  as  constant  values  of  p .*  Fig- 
ure 1 shows  the  details  of  the  transformation. 

8.  It  can  he  shown’'  that  a conformal  transformation  satisfying 

the  above  constraints  is 

N 

x(5,n)  = C + (B  sinh  nkn  + C cosh  nkn)  sin  nkC  (l) 


yCCjh)  = B + n + / (B  cosh  nkn  + C sinh  nkn)  cos  nk?  (2) 

o / j n n 

n=l 

w'here 

k = tt/X 

0 < X < X •,  0 < 5 < X (3) 

- B £ n £ e 

The  N values  of  B and  C and  B constitute  the  hi-curve  fitting 

conformal  mapping  coefficients.  These  coefficients  are  determined  by 

matching  the  prototype  coastline  and  sea  boundary  curves  at  p = B , 

respectively,  B also  being  a parameter  to  be  determined.  A.n  iterative 

3 

procedure  is  required  for  determining  B and  the  coefficients. 
Essentially,  the  procedure  is  terminated  after  convergence  of  B and 
the  coefficients  or  after  the  transform-generated  coastline  and  sea 
boundary  curves  are  in  good  portrayal  of  those  specified  (personal  '\riew 
of  the  user).  Typical  standard  deviation  between  transform-generated 
and  prototype  (slightly  smoothed  version)  over  the  length  of  the  curve 
is  of  the  order  of  2000  ft.  This  is  achieved  with  lO  iterations. 


For  convenience,  symbols  and  umisual  abbreviations  are  listed  and  de- 
fined in  the  Notation  (Appendix  A). 
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[■ 

i; 

b.  IMAGE  SPACE 

Figure  1.  Conformal  transformation  mapping  prototype  space  into 

image  space 

[ 


N = 100  and  at  a cost  of  9 min  of  CDC  6600  central  processing  time. 
Furthermore,  the  conformal  mapping  procedure  need  only  he  accomplished 
once  for  any  large  stretch  of  continental  shelf,  Ihe  computing  grids 
appropriate  to  areas  within  that  region  are  generated  depending  on  the 
nature  of  the  particular  study,  using  the  mapping  coefficients  previously 
determined.  The  computer  cost  in  generating  the  grid  is  inconsequential 
to  the  personnel  time  required  in  the  grid  makeup  or  assigning  the  cell- 
averaged  fluid  depth.  The  cost  in  generating  a grid  covering  a rela- 
tively small  extent  of  coast  must  he  weighed  against  a large  grid  of 
similar  high  resolution  over  a long  reach  of  coastline  which  costs  more 
per  prototype  hour  and  even  more  so  when  numerous  simulations  are 
required. 
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. PAET  III:  MATHEMATICAL  MODEL 


Stretched  Shelf  Coordinate  Svsten 


9.  Consider  the  transform-generated  coastline  and  sea  "boundary 
curves  shown  in  Figure  2.  The  orthogonal  curvilinear  mesh  associated 
with  the  (C,ri)  coordinates  in  prototype  space  is  designated  the  shelf 
coordinate  system.  The  display  of  that  system  in  prototj.'pe  space  is 
shown  in  Figure  2a  and  its  image  in  Figure  2'd.  Notice  that  only  a 
portion  of  the  entire  shelf  mapped  area  is  su'bsecuently  employed  for 
long-wave  computations.  Furthermore,  the  calculations  cannot  readily  be 
performed  in  the  image  space  because  the  cells  are  unevenly  spaced. 
Centered  computations  are  a prerequisite  for  finite  differencing  of 
partial  differential  equations.  To  provide  an  evenly  spaced  computing 
grid  and,  at  the  same  time,  to  obtain  the  desired  spatial  resolution 
with  the  fewest  possible  computational  points  require  a second  trans- 
formation. This  transformation  preserves  the  orthogonal  property  and 
allows  for  the  independent  stretching  of  5 and  of  n . The  grid  re- 
sulting from  the  second  transformation  is  the  evenly  spaced  computing 
grid.  Figure  2c.  The  coordinate  system  is  termed  the  stretched  shelf 
coordinate  system  (S,T). 

10.  The  stretched  shelf  coordinate  system  is  generated  by  in- 
dependently transforming  the  C and  h axes  in  the  following  manner: 

a.  Given  the  nearshore  region  of  principal  interest,  the  values  of 
C along  the  coastline  are  determined  which  will  produce  a con- 
stant relatively  fine  increment  of  coastline  arc  length,  S 
In  this  area  of  prime  interest,  the  line  BC  in  Figure  2a,  ^ 

the  constant  increment  of  arc  length  is  equal  to  AS  . How- 
ever, regions  AB  and  CD  show  that  for  the  same  AS  as  above, 
there  is  a relative  expansion  of  the  increment  of  the  coastline 
arc  length.  The  functional  relation  between  C and  S is 


S = S(Sp(e))  ' (U) 

where  the  expansion  of  S^  with  respect  to  S is  specified 
by  an  (arbitrary)  expression  of  the  form 


(5) 


+ 


(S) 


I 


b. 


where  A , B,  and  C are  constants  for  each  expansion 
region  (k).  For  the  example  in  Figui’e  2,  there  are  three  ex- 
pansion regions.  The  ensemble  function,  k = 1,  2,  and 
3 , is  determined  such  that  at  interface  points  the  function 
is  continuous  and  has  continuous  first  derivatives. 

Along  a particular  isoline  of  6 , , the  values  of  n are 

determined  which  will  yield  a constant  change  in  the  time, 

AT  , required  for  a long  wave  traveling  at  the  local  free 
wave  celerity  to  proceed  from  the  sea  boundary  (-3)  to  the 
coast  (3).  The  long-wave  travel  time,  T , is  calculated 


T = / (6) 

" bJ.) 

1 

J. 

where  is  the  distance  along  C'  > g is  the  accelei-ation  j 

due  to  gravity,  and  D is  the  local  depth  relative  to  mean  j 

sea  level  for  a standard  basin.  This  procedure  provides  a j 

fine  grid  spacing  in  the  nearshore  and  a coarse  grid  near  ] 

the  deep  sea  boundary.  The  relation  between  n and  T i 

is  given  by 


T = T(T^(S^(n)))  (7) 

The  incremental  values  of  T are  determined  from  Equation  6 
subject  to  the  (arbitrary)  expansion  relation  of  T(T  ) 
which  is  the  coimterpart  of  Equation  5-  Actually,  the  rela- 
tion T(T^)  is  a convenience  (seldom  used  in  applica- 
tions) which  permits  an  additional  degree  of  freedom  in  ad- 
justing the  relative  spacing  between  isolines  of  n . In 
most  applications,  T = T and  the  value  of  AT  is  that 
which  divides  the  total  long-wave  travel  time  by  an  integer 
number  of  lines  of  n • The  selection  of  AT  relative  to 
the  coarse  deep  sea  spacing  is  based  upon  a compromise  for 
providing  adequate  resolution  of  the  hurricane  winds  and  the 
deepwater  surge  with  a ninimTim  number  of  points.  At  first 
glance,  it  would  appear  that  the  T axis  for  hydrodynamic  long- 
wave calculations  is  time;  however,  this  is  not  the  case. 

Scale  factors  relating  the  transformation  of  n to  T o.re 
involved  in  the  equations  resulting  in  dimensions  of  i.engtli 
for  the  T independent  variable. 


11.  The  stretched  shelf  coordinate  system  provides  a grid  system 
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with  a finer  resolution  near  the  coast  than  at  the  deep  sea  boundary. 

The  expansion  curve,  5(3^),  stretches  the  alongshore  reach  of  the  grid 
w'hile  maintaining  a finer  grid  in  the  area  of  principal  interest.  In 
this  manner,  an  economy  is  achieved  in  terms  of  the  number  of  grid  points 
required  for  long-wave  simulation.  However,  because  the  preferred  ex- 
pansion curves  dictate  the  locations  (in  prototype  space)  of  the  fluid 
depths  required  for  the  long -wave  calculations,  the  depth  field  must  be 
redefined  for  different  combinations  of  stretching  functions. 


Governing  Equations 

12.  The  vertically  integrated  form  of  the  quasi-linear  long-wave 
equations  in  a Cartesian  system  is  well  known. ^ The  appropriate  forms 
of  these  equations  in  the  stretched  shelf  coordinate  system  are 

Transport  (momentura)  equations 


3Q 

~Ti  “ ^ Is'  ^ ’^3  ~ *^3 


(8) 


8G 


£2.  3. 


at  Fv  aT  "^t  ' '^T 


(9) 


and 

Continuity  equation 


M 

at 


+ 


1 3 


(FQg) 


7 k 


0 


(10) 


2 

where  Q is  the  volume  transport  per  unit  width  (units  of  length  /time), 
T is  the  wind  stress  divided  by  water  density,  a is  the  bottom 
frictional  resistance  stress  divided  by  water  density,  f is  the 
Coriolis  parameter,  g is  the  acceleration  due  to  gravity,  H is  the 
sea  surface  elevation  relative  to  mean  sea  level,  D is  the  total 
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instantaneous  depth  of  water  (=H  - D ), 

o 


D the  water  depth  relative 


to  mean  sea  level,  and  is  the  hydrostatic  elevation  of  the  sea 

surface  corresponding  to  the  atmospheric  pressure  anomaly.  The  in- 
dependent variables  are  time  (t)  and  the  coordinates  (S,T).  The  terms 
F , y and  v are  variable  scale  factors  associated  with  the  trans- 
formations. It  can  be  shown  that  F is  nondimensional  and  given  by 

1/2 

F = II  +1^1  I (11) 


rhe 


and  V scale  factors  are 


..  as 

U = E 

as  as 

p 


(12) 


and 


as  aT 

an  n n 

as  aT  aT 

n n 


(13) 


It  can  be  shown  that  y(6S)  and  v(oT)  have  the  dimensions  of  length. 
2 

Furthermore,  F yv  corresponds  to  the  Jacobian,  cf  the  transformaticn 
in  the  sense  that 


A = 


Ih 

p. 


2 

F yv 


dS  dT 


(lU) 


where  A is  the  area  of  a closed  region  in  prototype  space  whose 
corresponding  area  in  the  computing  space  is  R (see  Figure  2 for 
details).  This  relation  assures  that  all  area  enclosed  by  the  limits 
of  the  curvilinear  grid  in  prototype  space  is  accounted  for  in 
the  computing  space  by  use  of  the  scale  factors. 

13.  The  kinematic  wind-stress  components  and  x^  are  re- 

lated to  their  x,y  component  counterparts  (x  , x ) at  a given  point  in 

X y 

prototype  space  by 


T_  = T COS  0 + T sin  0 
S X y 


= -T  sin  6 + T cos  0 
T X y 


where 


0 = tan 


-1  / 3y/9g\ 

\3x/95/ 


The  relation  between  the  wind  stress  and  wind  speed  at  a reference 
anemometer  level  (usually  taken  near  the  water  surface)  is  talcen  as 


T = K , K = — C . 

10  pc 


where  is  air  density  and  p^  is  water  density,  C^,  is  a non- 

dimensional  drag  coefficient,  and  is  the  wind  speed  at  an  eleva- 

tion of  10  meters  above  the  water  surface.  The  value  cf  K is  taken  as 


1.1  X 10 


, if  <_  13.58  knots 


1.1  + 2.5  (1  - 13.58A;^q)^  10“^  , if  > 13.58  knots 


3,8  7 

This  form  for  K was  used  in  previous  studies.  Recent  compilation 

of  various  results  relating  to  the  wind  speed  at  10  meters  is  shown 

in  Figure  3.  As  reference,  the  dashed  line  in  Figure  3 shows  the 

relation  based  on  Eouation  19  with  p /p  = 0.0012  . In  the  middle 

a w 

to  low  wind  speed  range  (say,  less  than  60  knots).  Equation  I9  is 
representative  of  the  data.  It  is  this  range  of  wind  speeds  that  were 

used  by  Reid  and  Eodine'^  in  studies  of  historical  storm,  surges  in 

3 

Galveston  Bay.  VJanstrath  , on  the  other  hand,  in  studying  three  his- 
torical storm,  surges,  each  in  different  coastal  regions,  computed 
coastal  water  levels  in  good  agreement  with  observed  ccr.aitions  using 


15 


I 

o 


o 


I 


Figure  3-  Mean  values  of  the  neutral  drag  coefficient  as  a function 
of  low  and  high  wind  speeds  at  lO-m  height  (from  Garratt,  1977)  and 

that  used  in  SSUJtGS  III 
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Equation  19  and  estimates  of  the  wind  furnished  hy  the  Hydrometeorologi- 
cal Section  (Hydromet)  of  the  liational  V/eather  Service  (IWS),  National 
Cceanic  and  Atmospheric  Administration  (NOAA).  Those  maximum  winds 
ranged  from  95  to  120  knots.  However,  the  current  understanding  of 
the  surface  winds  from  hurricanes  obtained  by  aerial  reconnaissance, 

anemometer  records  from  oil  nlatforms,  and  sophisticated  numerical 

8 9 

marine  boundary  layer  models  is  that  the  surface  winds  at  10  meters 
are  not  as  severe  as  tho\ight.  Consequently,  with  liberal  estimates  for 
the  r.'.a>:imum  winds,  a conservative  drag  coefficient  at  high  wind  speeds 
is  necessary.  The  procedure  does  provide  good  surge  resiilts  as  evi- 
denced by  these  and  pre'/ious  studies.  Actually,  all  previously  reported 
open-coast  siirge  models  are  subject  to  recalibration  given  the  recent 
findings  of  hurricane  surface  winds. 

lU.  The  forms  of  the  seabed  frictional  resistance  terms  are 


and 


where 


a 


a 


S 


T 


KoQ 


.2 


^0, 


(20) 

(21) 

(22) 


and  is  a variable  nondimensional  drag  coefficient  that  depends  on 
the  seabed  condition  and  water  depth.  For  typical  seabed  conditions. 


K . F(0.02?)!  ,23, 

° (l.J^9)'^ 

where 

D*  = 0.25  , if  B < 0.5 

D*  = 1.0  , if  0.5  1 D < 1.5  (2I4) 

D*  = 2.0  , if  1.5  £ D < 2.5 

and  all  depths  D and  D*  are  in  feet. 
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Recursion  Eauations 


15-  The  numerical  analogs  of  Equations  8-10  are  based  on  centered 
difference  approximations  of  all  terms.  The  algorithm  treats  the  time 
dependency  explicitly  and  employs  com.puting  lattices  as  shown  in  Fig- 
ure The  recursion  equations  require  two  time  levels  to  complete  the 
cycle.  H values  as  well  as  h'  and  x are  comouted  on  the  lattice  at 

i3 

time  level  n . This  lattice  also  contains  the  permanent  storage  data 
of  B , 6 , X , and  Y w'here  X and  Y are  the  x and  y grid 
coordinates  of  H points.  For  a reference  index  the  transports 

Qg  and  are  computed  at  AS/2  and  AT/2  , respectively,  from 

the  H point  in  the  positive  axes  direction.  The  transports  are  com- 
puted on  the  lattice  at  time  level  n+1  . For  sm.all  time  steps.  At  , 
it  is  inconsequential  that  x is  computed  at  n and  applied  at  n+1  , 
Furthermore,  x is  computed  at  H point  locations  and  averaged  with  its 
neighbor  in  the  S and  T directions,  respectix^ely,  for  determination 

of  x^  and  x . These  approximations  for  x are  a result  of  efficient 
D i 

utilization  of  computer  time  and  memory/  considerations  and  can  be  sho’.+n 
through  numerical  experiments  to  be  accurate, 

l6.  Consider  that  the  transports  are  krlo■^^^x  at  time  level  n-1 
and  the  H field  is  known  at  time  level  n-2  . The  recursion  formula 
for  interior  grid  H points  is 

r:(l,J,n)  = H(l,J,n-2) 

2At[Fy(l,J)Qg(l,J,n-l)  - F^(l-l,J)Qg(l-l,J,n-l)] 


2At 

[F^(l,j)Q;j(l,J,n-l)  - 

-1)Q^(I,J-I,n-1)] 

M^ATVj,(.T) 

[fj^(i,J) 

2 

(25) 


where  I , J , and  n indices  express  the  S , T , and  time 
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r 


coordinates,  respectively.  At  is  the  numerical  time  step,  7,,  , , 

and  F..  are  the  point  values  of  F from  Equation  11  determined  at 

n 'b 

0^  , and  H points,  respectively,  is  the  point  value  of  p I’rom 

Equation  12  at  (typical  H grid  points  in  S direction),  v„  is 

J-  n 

the  point  value  of  v from  Equation  13  at  T_  (tj-pical  H grid  points 

• * . - «J  ... 

in  T direction),  and  is  the  map  factor  relating  protot;/T)e  length 

(feet)  to  x,y  units.  The  simulation  is  initiated  at  n = 1 with  all 
transports  and  previous  H field  equal  to  zero. 

IT-  Given  the  H field  just  computed  at  level  n and  the  previ- 
ous transports,  the  recursion  formulas  for  interior  grid  transports  are 


C^(l,J,n+l)  = 


Q,^(l,J,n-l)  2AtfQ^ 


2AtgDg[H(l+l,J,n)  - H(l,J,n)  - :-:2(l+l,J,n)  + H2(l,J,n)] 


M^ASGgPu(l)7u;i,J) 


At 


j^Tg(l,J,n)  + T^(l+l,J,n: 


(26) 


and 


Q„,(l,J,n+l)  = 


Q.^(l,J,n-l)  2AtfQ^ 


2AtgD^[H(l,J+l,n)  - H(l,J,n)  - :-:g(l,J+l,n)  + Hg(l,J,n) 

M^ATG^V^.(J;F^(I,J) 


where 


+ ^T^(l,J,n)  + T^(I,J+l,n)j 

= 0.5[H(l+l,J,n)  - D^(I+1,J)  + H(l,J,n)  - D^(l,j)J 


(2T) 


(2S) 


20 


(29) 


= 0.5  H(l,J+l,n)  - D^(I,.J+1)  + H(l,J,n)  - D^(l,J) 


(30) 


(31) 


(32) 


(33) 


, in  a similar  manner,  is  the  point  value  of  p from  Equation  12  at 

S_  + AS/2  (typical  Crid  points  in  the  S direction),  v is  the 

J.  i3  V 

point  value  of  v from  Equation  13  at  Tj  + AT/2  (typical  CL  prid 

^ H ^ 

points  in  the  T direction),  Xg  find  are  the  wind  stresses  in 

the  S , T directions,  respectively,  which  are  computed  at  H points, 

S T 

and  are  the  seabed  drag  coefficients  computed  from  Equation  23 

with  D in  Equation  2U  given  by  and  , respectively. 

k3  .1. 


Eour.dar*/  Condition 


l8.  Eor  closure  of  the  system  of  equations,  the  open  sea  and 
lateral  boimdaries  require  specification  appropriate  to  the  hurricane- 
induced  surge  computation.  I’.ar.y  different  continuum  type  boxmdai-y 
conditions  are  reported  in  the  literature.  In  principle,  it  should 


1 


i 


not  natter  v’aich  of  these  is  used.  If  the  corr.puting  grid  is  of  suf- 
ficiently large  extent,  the  solution  in  the  central  part  of  the  grid 
near  the  coast  should  he  insensitive  to  the  particular  one  used.  A 
popular  lateral  hoiindary  condition  which  is  also  used  in  SSUHGE  III 
is 

3Qo 

~=0  (3M 


The  sea  boundary  condition  is  that  the  water-surface  elevation  is  placed 
in  equilibrium  with  the  atmospheric  pressure  anomaly 


H(l,l,n)  = Hg(l,l,n) 


(35) 


-Along  the  lateral  and  sea  boundaries,  the  transports  are  computed  using 
modified  forms  of  Equations  30  and  31.  The  modification  is  to  compute 
an  average  value  using  only  the  two  adjacent  interior  points. 

19.  It  has  long  been  recognized  that  the  coast  boundary  condition 
plays  a major  role  in  the  nearshore  hydrodynamic  solution  of  the  govern- 
ing equations.  Most  open-coast  shelf  models  (whether  curvilinear, 
rectilinear,  finite  element,  implicit,  or  explicit)  treat  the  coast  as 
an  infinitely  high,  continuous  wall.  The  use  of  numerical  models  in 
simulating  circulation  and  water  level  conditions  is  made  difficult  be- 
cause of  application  of  such  a condition.  The  coast  boundary  condition 
developed  and  incorporated  into  SSURGE  III  is  the  Finite  Height  Barrier 
Coast  (FHBC).  The  FIIBC  provides  for  the  overtopping  of  low-lying  land 
and  bay  communication  with  the  open  sea. 

20.  The  FHBC  condition  is  to  allow  for  potential  ponding  areas 
landw'ard  of  the  shoreline.  Figure  5 shows  the  conceptual  design  of  the 
barrier  coast  and  adjacent  bay.  The  bay  is  sho-,m  to  be  of  several  incre- 
ments of  AS  in  length.  The  coast  flooding  routine  permits  a volume  of 
water  to  be  transported  across  the  nom.inal  shoreline  which  is  dependent 
upon  barrier  heights,  predicted  water  level  at  the  coast,  channel  en- 
trance characteristics,  and  the  water  level  in  the  ponding  area,  H^  . 

The  routine  is  applied  at  the  H time  levels  as  shown  in  Figure  •+. 
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OVERTOPPED  BARRIER  I 
SUBMERGED  , AND  ) 


EXPOSED 
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A(-  = CHANNEL  CROSS-SECTIONAL  AREA 
Cd  = CHANNEL  DISCHARGE  COEFF 
Cj  = SUBMERGED  BARRIER  COEFF 

Cp  = EXPOSED  BARRIER  COEFF 
OF  OVERTOPPING 

= BARRIER  ELEVATION 
L^  = COAST  LENGTH  AT  Zi, 

H^=  BAY  WATER  LEVEL 
H = COAST  WATER  LEVEL 


Figure  5>  Finite  height  barrier  coast  boundary 


21.  A prediction/correction  method  is  used  to  comijute  H(l,JM,n)  . 
The  philosophy  is  taken  that  the  coast  harrier  as  viewed  by  the  lonj?  wave 
is  essentially  a wall  over  much  of  the  tine  of  surge  development,  re- 
flecting the  majority  of  the  long-wave  energy.  This  results  from  the 
consideration  that  long-wave  reflection  occurs  when  the  wave  experiences 
a significant  change  in  fluid  depth  over  distances  that  are  small  rela- 
tive to  the  wavelength. The  typical  nearshore  bathymetry  combined 
with  coastal  barrier  elevations  assures  that  the  coast  is  essentially 

a wall  which  leaks  water. 

22.  The  prediction/correction  method  is  to  first  assume  total 
reflection  at  the  boundary  and  then  to  correct  that  prediction  of  H 
based  on  a finite  height  coast.  In  this  manner,  the  method  provides  a 
smooth  transition  from  the  time-dependent  circ'urstances  where  the 
boundary  is  not  flooded  and  no  corrections  are  required  to  the  cata- 
strophic flooding  coast.  Other  relations  for  predicting  H at  the 
coast  than  that  of  total  reflection  were  considered.  In  particular, 

a water-surface  slope  projection  method  was  tested.  Results  of  these 
tests  show  the  method  to  be  inappropriate  for  surge  simulations.^  Other 
tests  imply  that  the  slope  projection  method  is  applicable  to  tidal 
(free  wave)  sim.iilations. 

23.  At  time  level  n , the  water  level  at  the  coast,  , is 

predicted.  If  the  water  levels,  and  , exceed  the  minimum 

height  of  the  coast,  then  a volume  of  water  is  transported  across  the 
nominal  coast  for  the  time  interval,  n-1  to  n+1  , 


Submerged  barrier 

''a  ' 


where 


"h  = - »b 


(35) 


(36) 

(37) 


2k 


Cg  is  a nondimensional  submerged  barrier  coefficient  taken  as  0J4  , 

k k+1  k k+1 

and  ' are  lengths  of  the  coast  at  elevations  7.^  and  , 

respectively,  and  the  direction  of  the  flow  is  determined  by  the  sign 
of  (.positive  weans  water  is  removed  from  the  coastal  cell  and 

placed  in  the  ponding  area).  Barrier  heights  per  AS  section  centered 
on  H points  are  discretized  into  unit  elevations  (k).  The  length  of 
the  coast  at  each  elevation  is  obtained  from  topographic  maps,  beach 
surveys,  etc.  The  total  vol’jme  of  water  crossing  the  submerged  coast 
is  the  summation  of  Equation  35  for  k = 0 , 1...K  where  K is  the 
•integer  of  the  lesser  of  H or  . C- >nsidering  large  ponding  areas 
and  high  minimum  coastal  heights,  this  volume  is  generally  small 
compared  with  that  which  is  transported  through  channels  or  which  over- 
topped the  exposed  (on  one  side)  section  of  the  barrier.  These  submerged 
barrier  equations  reflect  the  assumption  that  frictional  effects  are 
dominant  in  the  bottom  layer  which  is  continuous  across  the  boundary. 

2U.  For  the  overtopped  harrier,  the  volume  of  water  crossing 
the  coast  is 

Exposed  (on  one  side)  barrier 


where 


''b  = 


h'^  - , if  < Z^  < h'*' 

b b b — 


(38) 


(39) 


Cp  is  a nondimensional  exposed  barrier  coefficient  taken  ns  0.2  and  the 
direction  of  flow  is  determined  by  the  sign  of  . The  total  volume 
of  water  overtopping  the  coast  is  the  summation  of  Equation  38  for 
k = 0,1,... M where  M is  the  integer  of  the  greater  of  or  . 


25.  The  vol\mie  of  water  associated  with  channel  communication  is 


V = + 2At  C g H - H,  + W D 
o — Dc  ' b'  cc 


SDc 


(»*0) 


where 


[h^  - R^l  , if  > 0 and  > 0 

, if  > 0 and  < 0 


, if  H <0  and  > 0 


(la) 


, if  11  <0  and 


< 0 


V/^  is  the  channel  width,  is  a nondimensional  channel  discharge 

coefficient  determined  for  each  entrance,  A is  the  channel  cross- 

c 

sectional  area  at  mean  sea  level,  and  the  sense  of  the  flow  is  taken 

toward  the  low  head  side.  An  estimate  of  C_A  is  possible  if  tidal 

D c 

observations  are  available  inside  and  outside  a bay  with  a constricted 
opening  to  the  sea 


C 6 = 

D c 


A 32-H 
s / o 


3gB^ 


(l42) 


where  A^  is  the  surface  area  of  the  bay  at  mean  sea  level,  T 


is 


the  tidal  period,  is  the  tidal  range  outside  the  bay,  and  is 

a nondimensional  bay  parameter  that  is  dependent  on  the  response 

or  phase  lag.  H.  is  the  tidal  range  inside  the  bay.  The  bay  parameter 

^ 12 
as  a function  of  response  or  phase  lag  is  shown  in  Figure  6. 

26.  Consider  that  values  along  the  coast  (uncorrected  for 

flooding)  are  known  at  time  level  n and  and  A^  are  known  at 

level  n-2  where  A^  is  the  bay  storage  area.  The  volume  of  water 
entering  (or  leaving)  the  ponding  s.rea  is  detennined  for  each  segment 
by  the  appropriate  sum  of  , and  . The  ratio  of  that 

volume  to  the  surface  area  of  the  grid  block  representative  of  , 

Aj  , provides  the  incremental  correction  to  the  predicted  coastal  water 
level  to  conserve  mass.  This  value  is  stored  in  H(l,JIi,n)  for  use  in 
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RESPONSE  FACTOR,  H/H 


NOTE:  H j = TIDAL  RANGE 

INSIDE  THE  BAY 

H = TIDAL  RANGE 

OUTSIDE  THE  BAY 

Figure  6.  Bay  parameter  Bp  as  a function  of  response  factor 

or  phase  lag 
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the  transport  computations  at  level  n+1.  The  nev  value  of  at 

time  level  n is  at  n-2  plus  the  incremental  vater-level  change 

from  the  instantaneous  and  even  distribution  of  the  entering  volume 
(the  sum  for  all  appropriate  AS  segments)  on  at  level  n-2  . 

From  II^  at  n , a new  bay  storage  area  is  obtained.  The  above  is 
repeated  for  each  ponding  area  and  its  designated  coastal  segments. 

Hiirricane  Wind  and  Atmospheric  Pressure  Models 


27.  The  hurricane  wind  and  atmospheric  pressure  models  as  given 
13 

by  Jelesnianski  are  employed  in  the  surge  studies  with  some  minor  ad- 
ditions. The  x,y-wind  components  for  a stationary/  storm  at  the  H grid 
points  are 

w = — [-(X  - X ) sin  (|)  - (Y  - Y ) ccs  (}.]F(r„)  (li3) 

X r^j  g e g e H 

and 


V = 

y 


- x^)  cos  <1)  - (Y^  - Y^)  sin  o]F(rj^) 


where 


rn  = [(X^  - x/ MY^ 


ihh) 


('i5) 


W is  the  stationary  storm  maximum  wind,  ({>  is  the  wind  ingress  angle 

reflecting  the  inward  flow  relative  to  that  wind  vector  tangent  to  the 

isovel,  P...  is  the  distance  from  the  storm  center  (X  ,Y  ) to  region  of 
’ H e e 

maximum  winds,  and  X and  Y are  the  (x,y)  coordinates  of  H grid 

S S . 
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points.  The  translation  of  the  storm  j^rovides  an  alteration  in  the  vind 
field.  The  x,y  supplemental  vind  components  due  to  storm  translation 
are 


= (T^  CCS  + Ty  sin  i;j)G(rjj) 


and 


vhere 


t = 


(-Tx  a in  + T^  cos  i>)g(iyJ 


- Cl  - (f)  - 90 


K_  , if  rjj  < Rj, 


(>+6) 


(hr) 


(h8) 


G(r,i)  = 


(it9) 


R,7Ti7  > 

n h 


T'x  and  T^  are  x,y  components  of  the  forvard  speed  of  the  storm  and 
n is  a rotation  angle  used,  primarily,  to  fine  tune  the  alignm.ent  of 
maximum  winds.  For  standard  operation  of  the  wind  model,  a is  set 
equal  to  90°  + (j)  . The  x,y-wind  stress  components  for  a m.oving  storm 
are 


= k(<4  * ^if) 


1/2 


w 


(50) 


and 


vnere 


(2  2\ 

= K I w + W W 

\ ^ yj  y 


/•  = w + t 

XXX 


W = V + t 

y y y 


(51) 


(52) 

(53) 
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and  K is  given  by  Equation  19-  The  stress  components  in  the  stretched 
shelf  coordinate  system  are  determined  by  applying  Equations  50  and  51 
in  Equations  15  and  16.  In  an  alternate  manner. 


and 


= W cos  0 + !■/  sin  6 
& X y 


W_  = -W  sin 
T X 


e + 


W cos 

y 


(510 


(55) 


where  6 is  given  by  Equation  17  and  and  are  the  wind  compo- 

nents in  the  (S,T)  system.  The  stress  components  are 


and 


(56) 


(57) 


28.  The  surface  atmospheric  pressure  field  associated  (but  not 
dynamically  coupled)  with  the  hurricane  wind  is 


P = P + (P 
o 


- P )e 
o 


(58) 


where  P^  is  the  central  press\ire  and  P^  is  the  far-field  pressure. 

It  is  often  observed  that  relatively  high  winds  remain  along  the  coast 
after  the  storm  has  proceeded  inland.  Through  numerical  experiments, 
it  is  found  that  reasonably  good  comparison  between  observed  (Hydromet) 
and  computed  winds  is  obtained  by  setting  to  be  the  distance  the 

storm  center  is  from  the  coast,  specifying  Wp  to  be  desired  w'inds  at 
the  coast  and  setting  a (measured  clockwise  from  the  storm  movement) 
such  that  the  maximum  wind  region  is  along  the  coast.  This  procedure 
allows  slightly  longer  simulations  to  be  performed  after  landfall  than 
would  be  permitted  otherwise  due  to  poor  wind-field  portrayal.  However, 
Rp  in  the  pressure  expression  (Equation  58)  is  not  allowed  to  increase. 
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If  this  procedure  is  followed  in  the  numerical  program  (it  is  optional), 
is  held  consiant  at  a value  RHIT  in  nautical  miles  (n.mi. ) after  the 
time  of  storm  landfall  (THIT)  vrhich  is  input  in  hours  after  the  start  of 
the  simulation. 

29.  The  various  program  parameters  that  are  used  to  generate  the 
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wind  and  atmospheric  pressure  fields  for  Hurricanes  Flossy,  Hilda, 
Betsy,^^  and  Carmen  are  given  in  Tables  1-U,  Only  for  the  postlandfe.ll 
winds  from  Hurricanes  Hilda,  Betsy,  and  Carmen  were  the  normal  routines 
modified  in  the  above  manner.  This  is  seen  by  noting  the  number  for 
THIT.  For  comparison,  the  computed  and  observed  (Hydrcmet)  i/ind  field 
at  selected  times  are  presented  in  Plates  1-3  for  Hurricane  Hilda  and  in 
Plates  for  Hurricane  Betsy. 

30.  The  terra  in  Ecpiations  8,  9,  26,  27,  and  35  is  computed 

in  feet  of  water  by 

Hg  = 0.0328(P^  - P)  (59) 

where  P , named  PINF  in  Tables  1-U,  and  P are  in  pressures  of 
millibars . 

31.  For  historical  storm.s,  all  necessai-y  input  parameters  can  be 

determined  if  Hydromet  has  assessed  the  maximum  winds  following 

the  same  procedures.  This  procedure  probably  results  in  a liberal  esti- 
mate for  VL  . Storm,  surge  simulations  for  forecasting  nurposes  typi- 
cally  involve  storm,  parameters  of  track,  P^  , P^  , and  R^^  . The 

only  other  necessary  input  for  this  wind  model  is  W . This  can  be 

n 

determined  from  the  correlation  of  a few  severe  historical  storms  w'here 

Hydromet  has  estimated  and  that  value  which  is  predicted  by  the 

^ 17 

Standai'd  Project  Hurricane  (SPH)  method  using  only  the  required  fore- 
cast storm  parameters.  A cursorj’-  examination  shows  that  the  ratio  V/^ 
to  V/^^^(SPH)  is  probably  in  the  range  1.2  to  1.3.  This  approach  is 
suggested  only  as  a temporary  solution.  Actually,  a dynamic  marine 
boundary  layer  hurricane  model  is  needed  that  is  independently  cali- 
brated to  several  storms  of  record  and  verified  against  even  more 
hir.torical  storms.  This  wind  model  is  comparable  to  the  treatment  of 
the  jihysics  of  the  atmosphere  as  SSURGE  III  is  to  the  hydrodynamics. 
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PART  IV:  HURRICANES  FLOSSY,  HILDA,  BETSY,  AI.’D  CARI-IEN 

SURGE  STUDY 


3B.  Hurricanes  Flossy,  Hilda,  Betsy,  and  Camen  affected  the 
Louisiana  coast  from  Atchafalaya  Bay  to  the  Mississippi  River. 

Hurricane  Flossy  is  the  only  storm  in  this  study  that  paralleled  this 
coastal  region  with  subsequent  landfall  near  Pensacola,  Fla.  Plate  T 
shows  the  location  nap,  the  four  storm  tracks,  and  the  computing  grid. 

The  Mississippi  River  delta  region  and  the  adjacent  levees  are  simulated 
as  nonovertopping  barriers  protruding  from  the  coast  in  the  manner  shown 
in  Plate  7-  For  three  of  these  hui’ricanes,  Kydromet  has  estimated 
the  surface  wind  fields  at  various  storm  locations.^  ,15»1<J  Information 
about  Hurricane  Carmen  is  from  personal  communication  with  Dr.  Joe 
Pellissier,  National  Hurricane  Center  (NHC),  NOAA,  Miami,  Flo.  It  is 
noted  that  the  storm  tracks  as  observed  from  Hydromet,  other  weather 
bureaii  offices  as  noted  in  poststorm  damage  survey  reports  of  the  U.  S. 

Army  Engineer  District,  New  Orleans  (NOD),  and  from  the  NHC  in-house 
historical  storm  track  record  are,  in  most  cases,  som.ewhat  different — 
in  particular,  the  track  of  Hurricane  Flossy  as  it  crossed  the 
Mississippi  River  delta  and  Hurricane  Carmen  as  it  approached  Lake 
Cliarles,  La.,  in  which  the  storm  was  decaying  (filling)  rapidly.  The 
stoim  tracks  presented  in  this  report  are  from  the  NHC  record  except 
that  the  track  of  Hurricane  Hilda  prior  to  landfall  is  moved  6 n.mi. 
to  the  east  and  the  track  of  Hurricane  Flossy  in  the  Mississippi 
River  region  is  a compromise  between  the  NOD  report  and  NHC. 

33.  Six  discrete  ponding  areas  with  storage  area  information  and 
FHBC  data  are  presented  in  Table  5-  Flood  Region  I is  from  I'arsh  Is. 
to  Point  Au  Fer  Is.  and  boiinded  inland  by  elevated  railroads  or  highways, 
spoil  banks  of  canals,  bayous,  or  other  waterways.  Flood  Region  II 
is  from  the  tei-minal  point  of  Flood  Region  I to  ’.line  Is.  and  bounded 
inland  by  the  spoil  banks  of  the  Homa  Navigation  Canal  to  the  east  and 
Bayou  Chene  to  the  west  and  highway  U.S.  90  to  the  north.  Flood 
Region  III  is  from  Wine  Is.  to  Belle  Pass  and  bounded  inland  by  the 
Houma  Navigation  Canal,  Bayou  Lafourche  to  the  east,  and  Bayou  Blue  and 
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a portion  of  U.S.  90  to  the  norch.  Flood  Region  IV  is  from  Belle  Pass 
to  Barataria  Pass  and  hounded  inland  by  Bayou  Lafourche,  Barataria  Bay 
Waterway,  and  U.S.  90  to  the  north.  Flood  Region  V is  from  Barataria 
Pass  to  the  vest  bank  of  the  Mississippi  River  levee  and  bounded  inland 
by  the  Barataria  Bay  Waterway  which  closes  with  the  levee  near  New 
Orleans,  La.  Flood  Region  VI  is  from  the  east  bank  of  the  Mississippi 
River  levee  to  the  entrance  of  Lake  Eorgne  following  the  seaward  o'Xtent 
of  the  marsh  area  and  bounded  inland  by  New  Orleans  and  U.S.  90.  All 
the  above-mentioned  bayous,  canals,  etc.,  have  spoil  or  natural  alluvial 
banks  associated  with  them  that  would  have  the  tendency  to  limit  flows 
from  one  flooa  region  to  another. 

3^.  For  Hurricane  Flossy  which  occurred  in  September  1996,  the 
coastal  surge  envelope  for  a wall  and  flooding  coast  boundary  condi- 
tions is  presented  in  Plate  8.  Phe  cor.parison  of  computed  water - 
level  envelope  and  observed  high-water  conditions  is  presented  in 
Table  6 for  coastal  areas,  '^‘able  7 shows  the  peak  value  of  and 

the  observed  conditions  for  each  flood  I'egion.  The  observed  conditions 
are  taken  from  the  NOD  poststor.m.  survey^^  which  is  presented  in  Plate  9. 
The  comparison  of  computed  and  observed  conditions  is  not  ns  good  as 
should  bo  expected. 

35.  The  term  "coastal  surge  envelope"  refers  to  the  peak  surge 
water  elevations  along  the  coast  without  regard  to  the  time  of  its  oc- 
currence. Vrhereas,  the  computed  coastal  water-level  envelope  is  the 
peak  water  level  from  the  time-i-ecord  r.ade  by  the  superposition  of  the 
expected  tide  and  the  surge  hydrographs  at  a particular  locatioii  wit'nout 
regard  to  the  time  of  its  occu''’rer.ce.  It  is  nore  difficult  to  compare 
computed  11^  ana  observed  irfiarii  high-vater  levels  because  the  tide 
cannot  be  simply  added  to  the  computed  sui-ge.  Furthermore,  loses 

some  meaning  in  surge  simulations  where  for  inlo.ni  locations,  the  local 
wind  setup  is  of  major  ii;portance.  Consider  a surge  of  5 ft  occurring 
at  a time  of  low  tide  which  is  expected  to  be  -5  ft  relative  to  msl. 
tine  surge  model  wo’.ild  compute  flc-cding  relative  to  barrier  elevations 
based  on  msl  and  surge  levels  ( rel  msl).  The  res-ult  is  an  erroneoiis 
computation  for  and,  to  a ’’ os.ser  degT’ec,  a.lso  for  the  coastal 
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surge.  However,  for  small  tide  ranges  such  as  exiifb  in  the  Gulf  of 
Mexico,  the  computations  for  the  open-coast  water  level  are  representa- 
tive of  the  superposition  of  tide  and  surge.  . This  prohlem  can  he,  in 
general,  avoided  if  the  tide  function  is  input  at  the  sea  boundary  and 
dynamically  computed  along  with  the  surge.  It  is  reconimended  that  this 
should  be  pursued  when  a generalized  predictive  tide  function  appropri- 
ate to  the  edge  of  he  continental  shelf  becomes  available. 

36.  Hurricane  Hilda  occurred  in  October  196^.  Plates  10-15 
show  the  water-surface  topography  together  with  the  depth- averaged 
’.'elocity  field  for  a wall  and  flooding  coast  boundary  in  6-hr  intervals 
over  a 30-hr  period.  The  role  of  the  I'ississippi  River''delta  in  limit- 
ing the  free  transport  in  the  continental  shelf  waters  is  clearly  seen 
in  the  series  of  the  above  snapshots.  The  effect  of  a flooding  boundary 
relative  to  a wall  is  seen  to  enhance  the  onshore-directed  transport  of 
water  in  the  nearshore  and  not-so-nearshore  region.  This  would  result 
in  more  sea  sediment  and  suspended  material  transported?  to  (and  over) 
the  coast  than  that  computed  in  ncnflooding  open-coast  models.  The 
principal  surge  area  occurs  at  Caillou  Bay  as  seen  in  Plates  12  and  16, 
which  show  the  surge  envelopes.  The  coastal  barrier  in  this  area  is 
typically  a beach  berm,  generally  less  than  3 ft  in  height  (m.sl).  Con- 
sequently, significant  flooding  occurred  from  I8OOZ  (Z  referring  to 
Greenwich  Mean  Time),  3 October,  to  approximately  0300Z  the  next  day. 

The  winds  at  0300Z,  h October,  are  directed  more  alongshore,  forcing 
the  inland  surge  to  the  east.  The  highest  observed  w'ater  level  of 
9.8  ft  occurred  inland,  to  the  east  of  Caillou  Bay  along  the  Houma 
navigation  Canal.  It  is  seen  that  the  difference  in  surge  envelopes 
for  a wall  and  flooding  coast  conditions  (Plate  16)  is  at  most  2 to 
3 ft.  This  difference,  although  seemingly  small,  reflects  a significant 

i 

I'looding  condition.  In  proper  perspective,  the  volume  of  water  flood-  j 

ing  the  coast  is  inconsequential  to  that  volune  entrained  in  the  near- 
shore open-sea  circulation.  Consequently,  water  can  be  removed  from 
the  coast  and  be  almost  instanta.neously  replaced  vfithout  significant 
effect  on  the  nearshore  water  level.  The  circula-tion  pattern,  however, 
is  greatly  changed  due  to  the  flooding. 
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37*  Tables  8 and  9 compare  the  computed  and  observed  high-vater 
levels  for  coastal  and  inland  regions,  respectively.  The  observed 
water-level  record  and  that  computed  (superposition  of  surge  and  ex- 
pected tide  as  determined  from  the  National  Ocean  Survey  tide  table) 
at  various  coastal  locations  are  presented  in  Plates  17-19*  An  oscilla- 
tion of  8 hr  in  period  is  shown  in  the  computed  hydrograph  at  Biloxi, 

Miss.  (Plate  19)*  The  pieriod  of  oscillation  is  characteristic  of  an 
edge  wave.  Multioscillations  occur  from  the  trapping  of  energy  between 
the  right  lateral  boundary  and  the  protruding  Mississippi  River  levee. 

However,  other  computed  water  levels  in  this  region  do  not  exhibit  such 
oscillations  ns  shown  in  Plates  l8  and  20.  Plate  17  shows  another 
typical  problem  in  open-coast  surge  computations — that  of  proper  wind 
specification.  Before  storm  landfall,  the  analj-'tically  computed  winds 
at  Eugene  Is.  are  always  directed  offshore,  resulting  in  the  enhanced 
drawdown  of  the  water.  Moreover,  the  islo.nd  is  on  the  leeward  side  of 
the  coast  and  land  frictional  influence  on  the  wind  is  not  considered 
in  the  hurricane  model.  Further  complication  occurs  in  light  of  recent 
developments  which  shows  that  the  surface  winds  in  the  forwa'^d  section 
of  the  storm  may  be  oubflowing  (a  negative  ingress  angle). The  re- 
sults presented  in  the  comparison  tables  and  in  Plates  17-19  are  in 
good  agreement  with  the  observed;  however,  improvement  should  be  ex- 
pected with  a better  wind  model.  Plate  21  presents  the  observed  high- 
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water  levels  as  reported  in  the  poststorra  survey. 

38.  Hurricane  Betsy  occurred  in  September  19^5.  The  water- 
surface  topography  and  velocity  field  comparing  a wall  and  flooding 
coast  boi'.ndary  are  shown  in  Plates  22-26.  The  effects  of  the  flooding 
boundary  on  the  nearshore  water  levels  and  currents  are  most  draitiati- 
cally  demonstrated  in  Plate  23.  Without  overtopping,  the  wall  condition 
is  tre.pping  water  in  excess  of  l8  ft  at  the  corner  between  the  east 
Mississippi  levee  and  coastal  wall.  With  flooding  of  the  marsh  area, 
only  12  ft  of  water  is  computed  along  the  levee.  The  coastal  surge 
envelope  (Plate  27)  clearly  shows  the  effect  of  a flooding  boundary 
versus  that  with  a wall.  The  velocity  fields  are  considerably  different 
with,  tigain,  the  enhancement  of  the  onshore  velocity  for  the  flooding 
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boundary.  Plate  23  also  shows  that' the  surge  is  just  beginning  to 

’ I 

develop  on  the  western  side  of  the  Mississippi  River  levee.  Plate  2^, 

3 hr  later  than  that  snapshot  in  Plate  23 j shows  the  developed  surge  in 
this  area  while  to  the  east,  sustained  relatively  high  water  is  con- 
tinuing. Plates  25  and  26  show  the  decreasing  surge  and  the  reentrance 
of  the  water  from  the  ponding  areas. 

39.  Tables  10  and  11  compare  the  computed  and  observed  high-water 
levels  for  coastal  and  inland  regions,  respectively.  Furthermore,  ob- 
served and  computed  hydrographs  at  two  coastal  locations  are  presented 
in  Plates  28  and  29.  The  hydrograph  location  in  Plate  28  is  to  the  left 
of  landfall.  The  problems  associated  with  the  computed  large  drawdoim 
have  been  discussed  previously.  No  presurge  water-level  anomaly  is 
added  to  any  of  the  computed  results  in  any  study  reported  herein, 
although  after  viewing  the  observed  hydrograph  at  Biloxi,  it  may  seem 
appropriate.  The  results  from  this  surge  simulation,  more  than  in  any 
other  presented  in  this  report,  show  the  necessity  in  treating  the  coast 
as  a finite  harrier.  The  comparison  of  the  computed  a.nd  observed  water- 

level  conditions  0,re  in  excellent  accord.  Plate  30  presents  the  observed 
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high-water  levels  as  reported  in  the  poststorm  survey. 

ItO.  Hurricane  Carmen  occurred  in  September  197^*  Snapshots  of 
the  water-surface  topography,  together  with  the  depth-averaged  velocity 
field  for  a wall  and  flooding  coast,  are  shown  in  Plates  31-33-  The 
principal  surge  developed  between  Caillou  Bay  and  the  entrance  to 
Terrebonne  Bay,  as  shown  in  the  surge  envelope  (Plate  3^),  at  a.pproxi- 
mately  0900Z,  8 September  (between  the  time  of  Plates  31  and  32).  The 
surge  does  not  develop  in  Atchafalaya  Bay  until  l800Z.  After  this 
tim.e,  the  forward  motion  and  location  of  the  storm  are  in  question  due 
to  its  rapid  decay.  Since  the  surge  developed  in  East  Cote  Blanche 
Bay  during  the  decay  stage  of  the  storm,  it  is  viewed  as  fortuitous 
that  the  peak  com.puted  and  observed  water  levels  are  in  good  agreement 
(Plate  35).  The  phasing  of  the  peak  water  levels  at  Li.ike's  Landing, 

La.  , is  in  poor  agreem.ent  probably  due  to  the  storm  proceeding  faster 
on  the  track  than  that  used  in  the  simulation. 

^il.  Plates  36-38  present  computed  and  observed  hydrographs  at 
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other  coastal  locations.  The  observed  hydrographs  in  the  region  to 
the  east  of  the  Kississippi  River  (Plates  3T  and  38)  show  a consider- 
able initial  water-  ei.evation  at  the  start  of  conputations . This  eleva- 
tion is  not  totally  a presurge  condition  because  the  action  of  the 
’..Lnd  is  seen  to  reflect  the  observed  peak  levels  some  20  hr  later.  On 
the  other  hand,  this  initial  condition,  ospeciai.ly  at  the  Gulf  Outlet 
Canal,  is  not  duo  to  the  hurricane  since  the  v/inds  o.a  the'  continental 
shel"?  are  very  weak  at  the  start- of  computations;  that  is,  if  the  surge 
si’r.ulation  had  begun  at  1200Z,  6 Reptember  (2U  hr  la  ad-vance  of  the 
.•tarting  time  in  Plai.e  76),. the  mo-.lel  would  not  have  shown  over  4 ft  of 
w-ater  at  the  canal  (.-ntrance  at  1200Z,  7 Goptember.  Other  processes 
are  occ-arring,  in-esimahly  not  assignable  to  the  storm,  sui-go,  for  what- 
eve-r  caused  the  trnnsi ent  water-level  bump  at  1700Z,  7 Ceptoir.ber 
(I’l.ate  37),  is  probably  respon-jible  for  a similar,  -^^ell-correlated  bump, 
at  2200Z,  7 September,  at  Biloxi  (Plate  38). 

h2.  Tables  . 12  and  13  compare  the  computed  and  observed  high- 
water  levels  for  coastal  and  inland  regions,  respectively.  Tliere  are 
no  observations  along  the  coast  in  the  principial  surge  area  for  com- 
parison with  the  compxited  value.  In  other  places  and  subject  to  the 
above  discxxssion,  the  comparison  between  computed  and  observed  water 

levels  is  in  good  agreement.  Plate  39  jiresents  the  observed  high-water 
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levels  as  reported  in  the  po.ststorm  survey. 
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PART  V:  C0MCLU3I0:iS  AI.'D  RECOMMil'.'DATIO'IS 


U3.  The  development  of  c.  two-dimensional,  time-dependent,  c'jrvi- 
linear,  open-coast  storm  surge  model  is  presented.  In  particular,  the 
model  treats  the  coast  boundai-y  as  a finite  height  harrier  which  is 
broken  with  bay  entrances.  The  I'lcoding  coast  routine  in  its  treat- 
ment of  submerged  and  exposed  (or.  one  side)  barriers  is  considered 
conservative. 

1*4 . Results  of  four  stonr.  surge  simulations  are  presented  in 
which  observed  and  computed  conditions  are  cc.mpared.  Model  results 
compare  most  favorably  with  ohosa  observed.  Ti'.e  finite  lioight  barrier 
approach  does  indeed  improve  the  h7dro;Iynamic  .monel's  behavior  in  the 
nearshore  region  of  the  open  coast  but  does  not  prove  an  accurate  pre- 
sentation of  inland  flooding  levels  where  local  wind  setup  is  of  major 
importance.  Other  techniques  are  required  utilizing  tlie  results  from 
the  open-coast  model  for  accurate  modeling  of  iriland  flooding  processes 
From  currently  available  data,  a relatively  low  value  for  the  wind  drag 
coefficient  at  high  wind  speeds  is  used.  This  is  necessary  to  compen- 
sate for  the  liberal  estimate  for  nhe  maxiraui/i  winds.  In  hindcasting 
situations,  it  is  not  guaranteed  that  the  sarae  procedure  is  used  for 
estimating  the  maxiirrira  wind.  The  wind  model  employed  in  these  surge 
simulations  from  historical  storms  requires  input  of  the  maximum  wind. 
The  problem  can  be  resolved  by  using  a wind  model  based  on  the  pl^sics 
of  the  marine  boundary  layer  that  will  provide  a wind  field  as  a result 
of  dynamic  meteorological  computations.  It  is  recommended  that  such  a 
model  be  obtained  (as  they  do  exist)  and  incorporated  into  the  liydro- 
dynamic  flooding  coast  model.  The  resulting  composite  model,  Dierefore 
v.'ould  be  appropriate  for  both  forecasting  and  hindcasting  and  involve  a 
drag  coefficient  reflective  of  the  currently  available  data. 

1*5 • The  following  recorrmendations  are  secondary  to  that  of  incor 
porating  a dynamic  wind  model  in  S3IJRGE  III.  It  can  be  shown  that  the 
astronomical  tide  does  distort  the  predicted  coastal  and  inland  flood- 
water  levels  when  it  is  not  incorporated  into  the  dynjimic  oom.putations ; 
that  is,  equivalent  results  are  not  obtained  by  the  superposition  of 
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surge  and  tide  and  that  simulation  where  the  tide  is  insludei  in  fr.- 
sea  boundary  forcing  condition.  For  the  east  coast  of  the  Vnited  f'ates 
w’lere  there  are  some  locations  with  large  tidal  ranges,  the  ride  sr.-uld 
be  included  in  the-,  computation.  Tlie  results  of  the  four  s'crm  surge 
simulations  performed  for  this  study  are  considered  sufficiently  acj,.- 
inate  using  superposition  of  surge  and  tide  since  they  occurred  in  Ghlf 
of  Mexico  waters  where  the  tidal  range  is  small.  It  is  rectr~ ended  that 
a generalised  tide  prediction  function  which  is  appropriate  ro  the  edge 
of  the  continental  shelf  be  Incorporated  into  the  surge  computations. 

h6.  It  is  recognized  that  other  approaches,  equations  and' coeffi- 
cients are  reported  or  can  be  developed  for  the  treatment  of  a f loo  ling, 
finite  height  barrier  coast.  These  should  be  investigated  and  the  more 
pi’omising  approach  tested  against  the  one  reported  herein. 
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Prograr’  Tire  of  0 hrs  = 1200  GMT,  9 Sept  1965 


Program  Time  of  0 iirs  •-  1200  GKI,  7 Sept  1974. 
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Table  6 

Coraparlson  of  Computed  and  Observed  Coastal 
High  Water  Levels  for  Hurricane  Flossy 


Coastal 

Location 

Computed" 

High  Water 

Level  (ft) 

Observed 

High  Water 
Level  (ft) 

6. 1 n.mi.  East  of 

East  Timbalier  I. 
at  Bay  Champagne 

2.3 

3.7  (EST* **) 

Southwest  Pass  Ent. 

2.0 

3.7 

Quarantine  Bay 

9.6 

10.8,  11.9,  12.1 

Pointe  A La  Kache,  LA  to 

Point  Pleasant,  LA 

8.9 

10.3,  10.5  (EST) 
10.7,  11.0  (EST) 

Breton  I. 

6.0 

7 to  8 (EST) 

i 


* Adjusted  for  astronomical  tide. 

**  Furnished  by  various  oil  companies  and  by  the  Freeport  Sulphur  Co. 
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Comparison  o 

Table  7 

f Computed  and  Observed  Inland 

High  Water 

Levels  for 

Hurricane  Flossy 

Computed* 

Observed 

Inland 

High  Water 

High  Water 

Location 

Level  (ft) 

Level  (ft) 

Barataria  Bay  to 
Mississippi  River  Levee, 
Flood  Region  5 

0.4 

3.6**,  3.6**,  5.2**,  5.5**, 
6.9**,  7.3**,  8.0** 

Lake  Borgne, 

Flood  Region  6 

2.5 

4.0,  5.2,  6.3,  8.5,  6.2t,  6.4t, 
6.9t,  7.7t,  8. It,  8.6t,  9.4t, 

10. 9t 

I 


NOTE:  Flood  Regions  1-4  not  affected. 

* Not  adjusted  for  astronoaical  tide. 

**  Observations  were  located  along  landward  side  of  coastal 
barrier. 

7 Observations  were  located  along  soutliwest  side  of  Lake  liorgne. 


Table  8 

Comparison  of  Computed  and  Observed  Coastal 
Hlsh  Water  Levels  for  Hurricane  Hilda 


Coastal 

Location 

Computed* 

High  Water 

Level  (ft) 

Observed 
High  Water 
Level  (ft) 

Gulf  Ent. , Freshwater  Bayou 

1.5 

2.3 

East  Cote  Blanche  Bay  at 

Lukes  Landing, 

2.5 

4.4 

Atchafalaya  Bay  at  Eugene  I, 

(w/o 

3.8 

flooding  coast,  4.3) 

3.3 

Lake  Felto,  2.4  n.ml.  Landward 
from  Ent.  to  Terrebonne  Bay 

(w/o 

6.8 

flooding  coast,  7.2) 

7.4 

Grand  I. 

(■w/ o 

2 

flooding  coast,  3.2) 

4.0 

East  Side  Barataria  Bay  near 
Gulf  Ent. 

5.0 

5.5 

Southwest  Pass,  Mississippi 
Pviver 

1.5 

3.2,  3.5 

Quarantine  Bay  at  Ostrica,  lA 

(w/o 

4.0 

flooding  coast,  4.4) 

4.6 

Ent.  Mississippi  River  - Gulf 
Outlet  Canal 

(w/o 

4.0 

flooding  coast,  4.4) 

4.6 

Biloxi,  MS 

(w/o 

: 5.2 

flooding  coast,  5.2) 

4.6 

* Adjusted  for  astronomical  tide. 


L 


Table  9 

Comparison  of  Computed  and  Observed  Inland 
High  Water  Levels  for  Hurricane  Hilda 


Inland 

Location 

Computed* 

High  Water 
Level  (ft) 

Observed 

High  Water 

Level  (ft) 

Atchafalaya  Bay, 

Flood  Region  I 

0. 

3.8,  5.1 

Atchafalaya  R.iver 
to  Houma  Canal, 

Flood  Region  II 

3.8 

None 

Terrebonne.-Timbalier  Bay, 
Flood  Region  III 

3.4 

3.0,  3.8,  4.0,  4.7, 

5.4,  6.6,  6.9,  6.9, 

7.0,  7.4,  9.8  1 

Bayou  Lafourche 
to  Barataria  Bay, 

Flood  Region  IV 

0.5 

None  - 

1 

} 

i 

i 

Barataria  Bay  to 

Mississippi  River  Levee, 

Flood  Region  V 

1.6 

1 

4.0  ! 

i 

i 

Lake  Borgne, 

Flood  Region  VI 

3.3 

3.8,  4.5,  4.7,  5.0, 

5.2,  5.3  1 

I 


* Not  adjusted  for  astronomical  tide. 


Table  10 


Comparison  of  Computed  and  Observed  Coastal 
High  Water  Levels  for  Hurricane  Betsy 


Computed* 

Coastal  High  Water 

Location Level  (ft) 


Observed 
High  Water 
Level  (ft) 


East  Cote  Blanche  Bay  at 
Lukes  Landing, 


1 

i 


Grand  I. 


East  Side  Barataria  Bay 
near  Gulf  Ent. 


West  Side  J'lisiissippi 
Iliver  Levee;  Empire,  LjI, 
to  Venice,  LA 


• 4.5 


3.8 


7.2  8.8 

(w/o  flooding  7.4) 


8.5  7.6 

(w/o  flooding  coast,  9.8) 


6.2  to  12.2  5.7,  7.4,  7.7, 

(w/o  flooding  coast,  8.8,  9.2,  10.4 
6.2  to  16.1) 


East  Side  Mississippi 
River  Levee;  Polnte  A La 
Hache,  LA,  to  Ostrica 
(Brenton  Sound) , LA 


12.4  to  15.2  13.6,  13.7, 

(w/o  flooding  coast,  14.4,  14.5,  15.7 

12.8  to  23.1) 


Gulfport,  MS 

9.0  (Approx.) 

10.7 

Biloxi,  MS 

7.6 

8.6 

Table  n 


Comparison  of 

Computed  and  Observed 

Inland 

lIiRh  Water 

Levels  for  Hurricane 

Betsv 

Inland 

Location 

Computed* 

High  Water 

Level  (ft) 

Observed 

High  Water 

T.evel  (ft) 

Atc.hafalaya  Bay, 

Flood  Region  I 

0. 

2.9,  3.9 

Atchafalaya  River 
to  Houma  Canal, 

Flood  Region  11 

0,9 

3.0,  4.3 

(along  Houma  Canal) 

Terrcbonne-Timbalier  Bay, 
Flood  Region  III 

1.8 

2.8,  3.0,  4.3 

Bayou  Lafourche 
to  Barataria  Bay, 

Flood  Region  IV 

0.5 

3.4,  5.4 

Barataria  Bay  to 

Mississippi  River  Levee, 
Flood  Region  V 

2,1 

3.4,  5.0,  5.3 

Lake  Borgne, 

Flood  Region  VI 

7.1 

6,4,  8.3,  8.8,  9.1, 
9.3,  9.8,  10.1,  10.7, 
11.0,  11.7,  14.4 

I * Not  .idjusted  for  astrorionical  tide. 

t 

t 
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Table  12 


(.oTnparlson  of  Computed  and  Observed  Coastal 
High  Water  Levels  for  Hiir r icane  Carmen 


Coastal 

Location 

Computed* 

High  Hater 

Level  (ft) 

Observed 
High  Water 
Level  (ft) 

East  Cote  Blanche  Bay  at 

4.3** 

4.5 

Lukes  Landing,  lA 

Atchafalaya  Bay  at 

4.5** 

4.9 

Lower  Atchafalaya  River 

Atchafalaya  Bay  at 

3.8 

3.7 

Eugene  I. 

(w/o  flooding  coast,  4.0) 

Grand  I. 

4.7 

4.2 

(w/o  flooding  coast,  4.8) 

Ent.  Mississippi  River  - 

5.1 

5.7 

Gulf  Outlet  Canal 

(w/o  flooding  coast,  5.4) 

Biloxi,  MS 

4.0 

4.5 

(w/o  flooding  coast,  no  change) 

Caillou  Bay 

9.2 

None 

(w/o  flooding  coast,  11.6) 

Ent.  Terrebonne  Bay 

8.7 

None 

(w/o  flooding  coast,  11.1) 

* Adjusted  for  astronomical  tide. 

**  The  time  of  high  water  occurred  after  1200  CST  8 Sept  1974.  The 
water  level  is  approximate  because  the  structure  of  the  storm  after 
1200  CST  is  uncertain  due  to  its  rapid  decay.  Computations  ended 
at  1800  CST  8 Sept. 
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Table  13 

Compaclsou  of  Computed  and  Observed  Inland 
High  Wafer  Levels  for  Hurricane  Carneri 


Inland 

Location 

Computed-'’ 

High  Water 

Level  .(ft) 

Observed 
High  Water 
Level  (ft) 

Atchafalaya  Ray, 

Flood  Region  I 

2.S,  3.2,  3.8, 
4.0,  4.8 

Atchafalaya  River  to 

Houma  Canal, 

Flood  Region  11 

3.8 

3.3 

Terrebonne-Timbalier  Ray, 

Flood  Region  III 

A.  4 

2.9,  3.8,  4.6, 

5.9,  11.6 

Bayou  Lafourche  to 

Barataria  Bay, 

Flood  Region  IV 

0.5 

None 

Barataria  Ray  to 

Mississippi  River  Levee, 

Flood  Region  V 

1.4 

3.1 

Lake  Borgne, 

3.4 

5.0,  6.0 

Flood  Region  VI 


* Not  adjusted  for  astronomical  tide. 

**  Computed  bay  level  rising  at  end  of  computation.  See  footnote *  ** 
Table  18. 
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